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INTRODUCTION 


The observation of the motion and structure of the 
magnetopause has been the topic of many papers in recent years 
(Holzer, et. al, , 1966, Heppner , g/l . al. , 1967 , Sonnerup and 
Cahill, 1967, Hyde, 1967, Sonnerup and Cahill, 1968, Cummings 
and Coleman, 1968, Kaufmann and Konradi, 1968). The importance 
of this topic lies in the fact that a clear understanding of 
the magnetopause would enable one to estimate the validity of 
the so-called open or closed model of the magnetosphere 
(Dungey, 1961, Beard, 1964), to check the various theories 
relating to the instabilities of the boundary (Sen, 1965, 

Parker, 1967 a and b, L'erche, 1967 , Southwood, 1968, Alpers, 
1969, Karlson, 1970) and to understand the origin of the drag 
at the magnetosphere boundary. This drag or transfer of 
momentum through the boundary is evidenced by the existence 
of the magnetotail and by the direction of the plasma flow 
at the boundary (Freeman, et. al., 1968, Warren and Freeman, 
1970) but there is no present agreement as to the exact physical 
mechanism involved in the drag (Piddington, 1960, Dungey, 1961, 
Axford, 1964, Siscoe, 1966, Tsuda, 1967, Eviatar and Wolf, 1968, 
Coleman, 1970, Cassen and Szabo, 1970). In this paper we 
present some observations of the boundary motion and structure 
at times of large drag and indicate the relevance of our 
observation to the problems described. 
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It was reported in a previous paper (Aubry, Russell and 
Kivelson, 1970, hereafter referred to as ARK) that on March 27, 
1968 the inbound 0G0-5 satellite recorded an inward motion of 
the magnetopause of about 2 Re in two hours. In this sequence 
of events, a reversal of the vertical component of the inter- 
planetary magnetic field was immediately followe'd by an inward 
motion of the magnetopause caused by an erosion of the magneto- 
sphere; magnetic flux was carried into the tail and the growth 
phase of a substorm began about half an hour later. The 
velocity of the inbound OGO-5 satellite roughly matched the 
average inward velocity of the boundary during about two hours, 
which enabled us to study the motion and structure of this 
boundary for a relatively long period of time. Moreover, 

0G0-5 was transmitting at its highest data rate throughout 
this period-; providing 56 vector samples of the field per second 
(sps), which allowed us to analyse the fine structure of the 
boundary . 

In section 2, information about the experiment will be 
provided; in section 3 some basic points of theory will be 
recalled; section 4 and 5 are respectively devoted to the 
oscillations and structure of the boundary. 

It will be concluded that our observations do not fit any 
steady open or closed model of the magnetosphere, that the motion 
of the boundary is at least not inconsistent with the Kelvin 
Helmholtz instability hypothesis, and that the existence of an 
observed variable magnetic field component normal to the boundary 



3 


could account for the drag responsible for the erosion of the 
dayside magnetosphere. 

2. THE EXPERIMENT 

2 . 1 The Orbit 

Some general comments about the 0G0-5 orbit were already 
included in a previous paper (ARK). For convenience. Fig. 1 
showing the projection of the 0G0-5 orbit in the GSM equatorial 
plane at the time of our observation is reproduced from this 
paper. Between 1700 and 1916 UT the distance between OGO-5 
and this GSM equatorial plane varied from 1.7 to 0.04 Re and 
so our measurements were performed essentially in this equatorial 
plane . 

2 . 2 The Magnetic Field Measurements 

The UCLA fluxgate magnetometer was designed to provide 
an accurate triaxial vector measurement of the magnetic field 
from perigee at low altitudes to apogee in the interplanetary 
medium. Each of the three orthogonal sensors has a dynamic 
range of +6^,000 gammas and in low fields can resolve field 
changes of l/8th of a gamma. This was accomplished with a 
basic magnetometer that measures +16 gammas and a set of coils 
that provide fields to null out the field at the basic magneto- 
meter to within +16y. The currents for these coils are provided 
in 64 steps of 16y and 128 steps of 1024y. This is accomplished 
as follows: If a field of greater than +16y or less than -16y 

reaches the basic magnetometer f this field is reduced in steps of 
16y until the field at the basic magnetometer is within its 


4 


operating range. When all available 16y steps have been applied 
(64 possible) a field of 1024y is applied and 63 of the 16y steps 
are removed. This stepping procedure has. a cycle rate of 500 hz 
which is far above the magnetometer sampling rates. 

The measured field consists then of three quantities: the 

number of 1024y nulling fields applied, the number of 16y nulling 
fields and the output of the basic magnetometer from +16 to -16y 
digitized in 256 parts, each l/8y. The sum of these three 
quantities for each independent axis gives the measured vector 
field . 

The basic magnetometer is operated as a closed loop magne- 
tometer with a frequency response that is flat to 150 hz and then 
rolls off at 20 db per decade above 150 hz . The three possible 
0G0-5 telemetry rates, 1, 8 and 64 kilobits per second, correspond 
to Nyquist frequencies of .43, 3.5 and 27.8 hz for the instrument. 
Since meaningful wave studies can be performed only if no signals 
above the Nyquist frequency reach the telemetry system, the 
output of the basic magnetometer enters a bit rate dependent 
filter before being digitally sampled. This critically damped 
fourth order filter has 8 db attenuation at half the Nyquist 
frequency, 20 db attenuation at the Nyquist frequency and 40 db 
attenuation at twice the Nyquist frequency. 

The satellite can simultaneously transmit data to earth 
(real time data) and store data on the spacecraft on magnetic 
tape for later transmission (playback data) . These data can 
be sampled at different rates: playback data is always sampled 
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at 1 kilobit per second whereas the real time data has three 
possible rates. Thus, the instrument has actually two outputs, 
each with its independent filter depending on the sampling 
rate of the digitization unit to which the signal is routed. 

The absolute accuracy of the measured field depends on 
many factors: the sensitivity of the magnetometer, the size 
of spacecraft fields and the possibility of drifts in the zero 
levels of the magnetometer. This magnetometer is the most 
sensitive fluxgate magnetometer ever flown on a spacecraft 
and is separated from the main body of the spacecraft by a 
twenty foot boom. However, there are other nearby experiments 
and this was a newly designed magnetometer. Comparing with 
the Goddard magnetometers on board which are on a similar 
boom restricted only to magnetometers, it was found that there 
was a slow drift from orbit to orbit of the apparent zero levels 
of the UCLA magnetometer. The Goddard magnetometers consist of 
a Rb vapor magnetometer and a fluxgate magnetometer both of 
which have been flown before on 0G0-1 and 0G0-3 and which 
provide consistent fields when compared. Data from these 
magnetometers, graciously supplied by the experimenters have 
been used to determine the zero level for each 0G0-5 orbit. 
However, in view of the slow drift of the zero levels the 
absolute value of any one component may be in error from one 
to two gammas. On the other hand the rate of drift is 
exceedingly slow compared to the time scales of physical 
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processes such as waves, discontinuities, etc. Thus changes in 
the field components can be accurately measured to the digiti- 
zation window of the experiment, 1/8 gamma. 

2 . 3 Particle Flux Measurements 

The pitch angle electron spectrometer was designed to 
measure the differential unidirectional flux of electrons in 
the energy range from 45 to 1200 kev . It consists of six 
approximately coplanar independent detectors, a pair aligned 
back to back and the other four spaced out through 90° relative 
to the back to back detectors. The detectors, each of which 
contains a plastic scintillator and a photomultiplier tube, 
make simultaneous measurements of electron flux in a cone of 
half angle 6°, and are capable of analyzing particle energies 
in eight logarithmically distributed intervals. The measure- 
ments are transmitted to earth (real time data) whenever 
possible, or alternatively, stored on tape for later trans- 
mission (playback data) . The maximum usable count rate corres- 

8 2 

ponds to a unidirectional flux of 10 electrons/cm /sec/ster. 

During the period of interest on March 28, 1968, the 
satellite telemetry rate was 64 kilobits per second, at which 
rate measurements are made in all detectors every .02 sec. 

The experiment was switched at ten minute intervals from a 
mode, the energy mode, in which all detectors repeatedly cycled 
through eight energy channels, to a second mode, the pitch angle 
mode, in which only total flux, in the energy range was measured. 
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The total flux in energy mode is obtained by averaging 

the differential flux over an integral number of cycles. As 

the averages used in this plot are taken over 70 cycles , the 

problem of short time variations is not critical. The statistical 

uncertainty (calculated on the basis of Poisson statistics and 

the digitization effects resulting from use of a prescaler) is 
3 2 

less than 10 electrons/cm /sec/ster in pitch angle mode and 
less than three times this in energy mode. 

3. MHD DISCONTINUITIES 

It is generally agreed that the magnetopause must be mainly 
a MHD tangential or a rotational discontinuity, the characteristics 
of which are shown in Fig. 2 (Landau and Lifshitz, 1960). In 
this figure, p, v and p are the density, flow velocity and 
pressure of the plasma; _H_ is the magnetic field, and subscripts 
n and t refer to the components normal and tangential to the 
boundary. The square brackets [ ] represent the difference 

between the quantity on the two sides of the discontinuity. 

The theory does not define the variation of the field 
inside the discontinuity. We know only that for a rotational 
discontinuity, and H 2 , the projections in the plane of the 
boundary of the magnetic field on either side, should be equal 
in magnitude, a normal component of H can exist, and the direction 
H^-H 2 = (H t 3 must be parallel to the vector [v fc ] . This last 
requirement is in fact a very useful test; indeed if we consider 
the equatorial region of the magnetopause, the plasma flow 
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velocity on either side must be mainly horizontal and so must ■ 
the vector [v ] . This implies that if the boundary corresponds 
to a rotational discontinuity the vector must be mainly 

horizontal. We shall use these relationships in -sections 5 and 6 
to identify the nature of the field discontinuity at the 
magnetopause . 

4. MOTION OF THE BOUNDARY 

The average inward motion of the boundary between 1700 and 
1900 has already been reported (ARK) . In that paper it was 
shown that boundary oscillations superposed on the main inward 
motion resulted in multiple crossings. We are interested here 
in studying these oscillations of the boundary in detail. 

First, Figs. 3 and 4 present our data. Figure 3, reprinted 
from ARK, shows the one minute averages of the magnetic field 
in the geocentric solar magnetospheric (GSM) coordinate system 
from 1500-2100 UT . The sequence of multiple crossings occurs 
between 1700 and 1916. The magnetic field and particle data 
are shown in Fig. 4 which consists of three panels. Each panel 
refers to the same time interval and contains forty minutes of 
data. The top line in each panel represents 4.6 second averages 
of B^, the component of the magnetic field along the axis 1 of 
the satellite reference system (in which reference system, the 
axes are labelled 1, 2, 3). For the period considered, this 
axis is nearly antiparallel to the Z GgM axis (the angle varies 
between 170° and 180°) and crossings of the magnetopause appear 
as reversals of which is positive in the magnetosheath and 
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negative in the magnetosphere. In order to avoid any ambiguity, 
these two regions are labelled in Fig. 4 for the first crossings. 
The middle line represents 10 sec averages of the particle fluxes 
for electrons with energy greater than 50 keV. These fluxes were 
computed from the readings of detector 5 which was oriented at 
the largest pitch angle. Flux greater than 4 x 10 particles/cm / 
sec/ster is indicated by a straight line, while gaps in the curve 
occur where data is not available. The bottom line of each 
panel represents the total magnetic field amplitude. 

The first clear encounter of the boundary took place at 
1700 UT (Fig. 3) when the geocentric distance of the satellite 
was 12.8 Re, (point A in Fig. 1). The field amplitude was about 
the same on both sides of the boundary; only the horizontal 
component varied. While no particles were detected before 1700, 
a detectable flux of energetic electrons appeared at 1700, and 
was observed till 1731 (see Fig. 4). This verifies that the 
satellite entered the magnetosphere at this time. This fact is 
confirmed also by the JPL solar wind experiment and the Lockheed 
ion mass spectrometer on board 0G0-5 . At this time the field 
in the magnetosheath was northward and it appears from Fig. 3 

i 

that there was essentially one single crossing. By looking at 
4.6 sec average data we checked that the small oscillations just 
before 1700 in Fig. 3 were in fact three partial crossings. So 
at this time the boundary was not perfectly steady; small 
amplitude oscillations of 1 to 3 minute periods were. present , but 
their amplitude was much smaller than that of the oscillations 
observed after 1730 UT. 
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At 1730 after the reversal of the field in the magnetosheath 
(as already discussed in ARK), the .sequence of multiple crossings 
began (Figs. 3 and 4). In order to interpret these crossings, 
we will refer to two models of boundary as sketched in Fig. 5. 
Model "a" presents a superposition of large amplitude oscillations 
with periods (T^) of several minutes, and small amplitude 
oscillations with periods (T^) of less than one minute. In model 
"b" the short and long period oscillations are again superposed 
but with approximately the same amplitude. 

From 1730 to 1840 UT (points B and G in Fig. 4), the 
occurrences of short duration multiple crossings appear at well 
spaced intervals and the observations can be interpreted in terms 
of model "a". The long periods, T^, range from 3.5 to 6 minutes 
and the short ones, T ^ > can be as short as 10 seconds. For 
instance, at 1807 a sequence of partial crossings occurs (the one 
labelled "2" in Fig. 4 will be studied in detail in section 5 
and will be shown to correspond to a 10 sec oscillation of the 
boundary). Between about 1840 and 1903, as discussed in ARK, 
the magnetosheath field turned more southward; this produced 
a new inward motion of the boundary and the multiple crossings 
observed between 1840 and 1906 present a quite different aspect 
from the former ones. Since now the short duration multiple 
crossings occur almost continually, they can be explained in 
terms of model "b" (Fig. 5), with a long period, T^, of about 
5 minutes and a short period, T^, 


of about 1 minute or less. 
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Some minutes after 1900 the magnetosheath field became less 
southward again and the two last crossings at 1912 and 1915 
reflect this situation. Namely the relative amplitude of the 
short period oscillation has somewhat decreased. 

Prior to the reversal of the magnetosheath field at 1730, 

the peak flux observed by the energetic electron spectrometer 

4 2 

is 4 x 10 particles/cm /sec/ster. After 1730 the peak flux 
measured inside the magnetopause increases by an order of 
magnitude. Oscillations in intensity of the electron flux with 
periods of several minutes, consistent with the model of large 
amplitude oscillations of the boundary occur in the regions 
inside the magnetopause between 1730 and 1844 and af.ter 1905. 

The period between 1844 and 1905 is characterized by rapid 
variations in the electron flux, frequently with periods of 
the order of one minute, and lower peak intensities. Although 
the oscillations of several minutes period as recorded by the 
particle experiment are in general agreement with the oscillations 
of the magnetic field, it can be checked from Fig. 4 that the 
relative position of the magnetic field and particles boundaries 
are variable: obviously the picture of a magnetic field 

boundary dividing two regions with- and without energetic particles 
is too simple. 

To be more specific about the motion of the boundary we 
must describe the motion of the wave forms sketched in Fig. 5. 
Namely, are they drifting from the front of the magnetopause 
toward the tail, are they standing, or are they coming from the 
tail, and what is their velocity? This problem is an important 
one, because to compute the thickness of the boundary from the 
time duration of the observed crossings, we must know the 
velocity of the boundary relative to the satellite. Quite 
generally, only the velocity normal to the quiet boundary is 
computed (Holzer, et. al., 1966, Smith and Davis, 1970). However, 
it appears from theoretical studies that the magnetopause 
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is often unstable relative to the Kelvin-Helmho ltz instability 
(Sen, 1965, Southwood, 1968) especially when, as in our example, 
the plasma flow velocity is normal to the magnetic field. It 
is, therefore, quite probable that the oscillations we observe 
are a consequence of a Kelvin-Helmholtz instability of the 
boundary. If this is so, it has been shown by Southwood (1968) 
that the motion of the boundary is produced by the superposition 
of two MHD waves propagating parallel to the boundary with 
opposite phase velocity on either side. If the tailward plasma 
flow velocity parallel to the boundary is U inside the magneto- 
sphere and U + V just outside the boundary, the motion of the 

V 

waveforms in Fig. 5 should be tailward with a velocity U + -j. 

In the next section we shall show that these multiple 
crossings can correspond to the tailward drift across the 
satellite of ripples such as sketched in figure 5 with a velocity 
of the same order as the plasma flow velocity. This is not 
inconsistent with the hypothesis that such oscillations are 
produced by the Kelvin-Helmholtz instability and indicates that 
the oscillations are not pulsations of the boundary as a whole. 

5. STRUCTURE OF THE BOUNDARY 

We shall now study some characteristic crossings in our 
two hour period of observation. For that purpose we shall use 
the eigenvalue technique discussed in Appendix 1. The magnetic 
field will he presented in the reference system xyz of the eigen- 
vectors: x (largest eigenvalue) is roughly vertical and is 
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arbitrarily oriented positive upward? z_ (smallest eigenvalue), 
normal to the boundary, is generally roughly horizontal and is 
arbitrarily oriented positive outward. As a consequence, 
(intermediate eigenvalue), roughly horizontal too, is positive 
tailward, and the xy plane is the plane of the boundary. 

The directions of the and x. axes corresponding to the 
smallest and largest eigenvalue for some of the crossings before 
1840 are shown in Fig. 6 in the GSM reference system; one 
second averages were used for computing these eigenvectors. 

Each point represents the tip of the projection of unitary x 
(circle) or z_ (arrow) vectors on the xy GSM plane. Circles 
corresponding to angles of 30°, 45° and 90° with the GSM Z axis 
are shown. The arrow labelled z and the circle labelled x 
correspond to the unperturbed boundary. The cross section of 
this boundary in the GSM equatorial plane has already been 
defined in ARK by the relation 

R = —5 

1 + e cosa 

R is the geocentric distance and a the sun- ear th- satellite 
angle (about 40°), e = 0.35. The satellite is not exactly in 
the GSM equatorial plane at the time of observation: in order 
to compute jc and _z, a magnetopause axisymmetric around the 
earth-sun line is assumed. 

As explained in Appendix 1, there can be a considerable 
scatter of the eigenvectors computed from different time 
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intervals even within one particular crossing. We caution the 
reader not to consider each individual direction in Fig. 6 as 
meaningful, but to examine the general trend of the pattern. 
Namely, the jc directions (circles) make angles generally less 
than 30° with the GSM Z axis; the _z directions are extremely 
variable roughly in the horizontal plane, so the waveforms of 
Fig. 5 would correspond roughly to an equatorial cross section 
of the boundary. This corresponds to a fluting of the magneto- 
pause and is in agreement with the results of Kaufman and 
Konradi (1968) . 

The sequence of crossings from 1700 to 1916 can be divided 
into three parts. The first one around 1700 when the magneto- 
sheath field is upward, is characterized by the absence of 
important fluctuations of the boundary (Fig. 3). The second 
interesting zone extends between B and C (Fig. 4) and represents 
a less stable state of the boundary following the partial 
reversal of the interplanetary field. The crossings numbered 1 
to 4 in this zone will be discussed more extensively. The third 
zone after C (Fig. 4) represents a more unstable state of the 
boundary following the total reversal of the interplanetary 
field as already discussed in ARK. The crossings labelled 5 and 
6 will be discussed. We shall try to determine whether or not 
the various states of unstability and the different drags 
observed in these three parts of the data are associated with 
a specific structure of the boundary. 
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5 . 1 First Zone 

In the first zone, the crossing at 1700 has been studied: 

when examined with the 1 sps data, it consists mainly of a 

rotation of the field lasting about 15 seconds with no normal 

-2 

component: the computed- average for B is less than 10 y and 

so less than the digitization error of the experiment. The 
crossing was then broken into portions of about 2 sec length 
and the eigenvalue analysis performed for each of them using the 
56 sps data. This analysis confirmed that the boundary had a 
rather steady orientation during the whole crossing. However, 
the components of the magnetic field normal to the boundary 
varied from 10y to -10y from one interval to the other. This 
could be due to the presence of waves inside the boundary. 

5 . 2 Second Zone 

Crossing #1 

The orientation of the eigenvectors for this boundary are 
shown in Fig. 6, where for convenience they are labelled x^ 
and z^. It appears that the z^ axis is very close to the 
direction of the unperturbed boundary normal. 

The variation of the three components of the magnetic field 
in the reference system of these eigenvectors is shown on Fig. 7 
with a sampling rate of 56 sps. The characteristic features of 
this crossing are: reversal of B x , maximum of B^ inside the 
boundary, compr essional oscillations of about 4 sec period and 
5y amplitude inside the boundary. The average component, as 
computed from the one second averaged data, is 2 . 4y . 
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The -variation of the tip of the one-second averaged 'magnetic field 
in the tplane of the boundary as seen from the earth "is ' sho^d Sii Fig. 8. 
This variation is mainly a rotation and the two dashed-line 

circles centered at C have been drawn for visual convenience. 

H H 

However the ratio 1 - 2 , which should be zero for a 



rotational discontinuity, is 0.3 and the direction is 

mainly vertical (see the direction of x^ in Fig. 6). These two 
facts refutes the hypothesis that crossing #1 is a rotational 
discontinuity, (cf. section 3). In order to check the 2 . 4y 
normal component, the crossing was then broken in intervals of 
2 sec length and each of them analyzed by using the 56 sps data. 
The computed normals to the boundary were rather scattered and 
the corresponding normal components ranged from +4 y • to — 23y. 
Even if we ignore this last value, there was so much scatter 
in the normal field that it is probable that the motion of the 
boundary forbids us from drawing any conclusion from the various 
normal components obtained and so we do not know if crossing #1 
corresponds to a tangential discontinuity. Our main purpose in 
presenting this crossing was to show an example of a rotating 
magnetic field across a boundary which is not a rotational 
discontinuity . 

Crossing #2 

This is a double crossing; the satellite which is first 
in the magnetosphere crosses through a bump in the boundary 
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and makes a short excursion into the magnetosheath lasting 
several seconds. To study this part of the data, we selected 
an interval of 21 sec which contains the spike labelled 2 in 
Fig. 4. Figure 9 shows the variation of the magnetic field in 
the reference system (1, 2, 3) of the satellite with a sampling 
rate of 56 sps , 

In order to study the magnetic field in a more convenient 
reference system, an eigenvalue computation from the one-second- 
averaged data was performed for this interval. This allows us 
to separate the B component which reverses at each crossing 
from the and components. We assume that V, the average 
direction of the field in the yjz plane is also the direction of 
the plasma flow which carries the frozen-in field, as well as 
various bumps on the boundary surface. Each second we plot along V_, 
starting in a, the direction of the field projection on ^z. (arrows) 
and the amplitude of the component (Fig. 10). Assuming that 
at any instant in time the amplitude and orientation of the 
field are constant along a direction perpendicular to V, the 
direction of the flow, we draw the projection of lines of force 
in the plane yz. The resultant sketch of the boundary explains 
simultaneously the variation of amplitude and orientation of the 
field. Between the points a and d the whole picture is self- 
consistent: namely, given the position of the first crossing 
at b and the orientation of the field during the following 10 
sec, one can predict the position of the second crossing at c; 
moreover, the values of the component of the magnetic field 
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around the crossing c can be predicted from their values 
around the crossing b. This is no longer true beyond the point 
d (Fig. 10) where a change in orientation and/or amplitude of 
the plasma flow probably occurs. 

For each crossing, b and c, the eigenvector and n c normal 
to the boundary were computed from the high time resolution 
data (56 sps). They are shown in Fig. 10 (bottom left). The 
orientation obtained for n^ is in agreement with the general 
shape of the boundary. On the other hand the agreement is not 
so good for n c . Owing to the great variability of the boundary 
at this time we did not try to obtain a better agreement for n c 
by selecting several time intervals. The variation of the 
magnetic field during crossing b, as seen in the reference system 
of the boundary is shown in Fig. 11. The whole crossing b lasts 
about one second between M and N. The reference system used is 
the one corresponding to eigenvectors for crossing b alone. 

Except for the large difference in time scale (1 sec instead of 
10 sec), the pattern is very similar to that observed in 
crossing #1 (Fig. 7) namely 

- reversal of B 

x 

- maximum of B inside the boundary 

y 

- fluctuations of B^ around some average value ('t-O.Sy) 

We think that the usefulness of the method we used for 

reconstructing the shape of the boundary is confirmed by the 
coherency of the results obtained. However, we must emphasize 
its drawbacks: The method would be exact for a rippled tangential 
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boundary embedded in plasma flow with a velocity uniform in 
space and constant in the interval of time of the chosen sample, 
here about 20 sec, and providing the frozen field hypothesis is 
satis'f ied . We know that all these conditions are not met; for 
instance, in the crossing c we obtain a normal component of about 
4y , and around the point d there is some change in the plasma 
flow. Finally the instantaneous direction of the plasma velocity, 
assumed to be parallel to V_, at about 30° from the unperturbed 
boundary for 20 sec, is justified only by the results obtained. 
Owing to these limitations. Fig. 10 can give only qualitative 
information about the shape of the boundary but we think that 
this information is useful. Namely, the very rapid oscillation 
of the field as exemplified by the spike labelled 2 in Fig. 4 
is consistent with the tailward drift across the satellite of 
small ripples on the boundary. 

It is possible to obtain some information about the distances 
involved in the model of Fig. 10. From the Explorer 35 measure- 
ments we know some parameters of the solar wind: temperature: 

T =9x10^ °K, velocity; V =470 km/sec, Alfven Mach number = 6.2. 
s s 

From the model of the plasma flow around the magnetosphere for 
an Alfven Mach number of 8 (Spreiter and Alksne, 1969) the 
value of these same parameters can be determined for the point 
of observation (equatorial plane, 0900 LT) ; we obtain 

V = 0.4 V = 200 km/sec 
s 

T = 19 T = 1.7xl0 6 °K 
s 
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Let us assume that during 13 sec, between the points a and d of 

Fig. 10, the velocity V of the plasma flux was oriented along 

0 V . This gives us a distance be, which is half a wavelength of 
the ripple, of 200 x 9 = 1800 km 'v 0.26 Re. What was the 
thickness of the boundary at this time? At b the satellite 

crosses the boundary in 1 sec, at an angle of about 45°. This 

gives a boundary thickness of 200 x 0.7 = 140 km. Now the 
Larmor radius for 1.7 x 10^ °K protons in a 40y magnetic field 
is 64 km. So the computed thickness of the boundary at b is 
about 2 Larmor radii. Theoretically this thickness is a 
minimum one for a reversal of the field and, consequently, the 
assumed velocity of the ripple is a minimum one too. This 
verifies that the ripple responsible for the spike #2 in Fig. 4 
has to move tailward with a velocity at least of the same 
order of magnitude as the plasma flow velocity. 

Crossing #3 

The position of this crossing in the whole sequence, is 
shown in Fig. 4. Crossings #3 and # 4 have been studied 
successively and in more detail than the former ones. Indeed 
they are less than 2 minutes apart which gives an opportunity 
to study the important changes occurring in the boundary during 
very short intervals of time. Moreover they seem to be "clean” 
with the time definition of Fig. 4 and this will allow us to 
show the importance of high sampling rates for studying the 
magnetopause. We wanted to check carefully on these crossings 
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the tailward drift hypothesis as well as the existence of 
variable components of the field normal to the boundary. This 
led us to perform several eigenvalue computations for each 
crossing; this is presented in the Appendix. 

In regard to the crossing #3, le-t us first examine the 
variation of the magnetic field across the boundary. We have 
to select some of the time intervals presented in the Appendix: 
first the interval n lx (Fig. A^ is chosen because it represents 
one second averaged data for the whole crossing (about 10 sec). 

In Fig. 12 the variations of ]5 in the reference system of the 
boundary (computed from the n^ interval data) is shown. The 
classical signature of the boundary appears: reversal of B^, 
maximum of B and oscillations in B around an average value 

y 2 

(5y) • These features are very similar to those of crossing if 1 
(Fig. 7) or 2 (Fig. 11). To look in more detail at the region 
of the field reversal, we choose the time interval (Fig. A^) 
which contains 56 sps data for this region. 

In Fig. 13 the variation of the projection of the magnetic 
field in the plane xy of the boundary (computed from the 
interval data) is shown as seen from the earth. ' This- 
Variation ’-should be compared with the rotation in Fig. 8 . ~ 

During the interval AB , the corresponding B z component varies 

-2 

between -2y and +2y with an average of 3.5 x 10 . If we compute 

the average B component over the whole interval, the result 
is somewhat different (-0.5y as indicated in Table 1 of the 

This result emphasizes the difficulty' of measuring 


Appendi-x) . 
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the component of the field normal to the boundary. It appears 
that the structure of this part of the crossing which lasts 
less than two seconds is rather complicated: a central section 
NP nearly force-free is squeezed between two intervals MN and 
PB where opposite Lorentz forces normal to the boundary are 
maximum (B^ and VB are nearly parallel in these two intervals) . 

As a matter of fact, the use of magnetohydrodynamic concepts 
should be considered here with caution because it will be shown 
that the total thickness of the boundary during crossing #3 
is of the order of 3 Larmor radii. 

In regard to the existence of a normal component of the 
magnetic field to the boundary, it is concluded in the Appendix 
that no clear result can be drawn for this particular crossing. 

Crossing #4 

The position of this crossing in the whole sequence is 

shown in Fig. 4; it occurs less than 2 minutes after crossing # 3 

and we will use for #4 exactly the same kind of analysis as for 

#3. A discussion of the eigenvalue analysis is presented in 

the Appendix. In order to study the variation of the magnetic 

field across the boundary we first select the interval ng which 

represents one second averaged data for the whole crossing 

(Fig. A^) • In Fig. 14 the variations of B^ in the reference 

system of the boundary (computed from the ng interval data) 

is shown. The reversal of B and the maximum of B are observed. 

x y 

We chose the interval N, (56 sps) to look in more detail at the 

b 
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region of the field reversal. In Fig. 15, the variation of the 
projection of the magnetic field in the plane of the boundary 
(computed from the interval data) is shown as seen from the 
earth; one observes mainly a rotation. During this interval 
the corresponding B component varies between 11 . 3y and 12. 6y 
with an average of about 12y. The variation of the field in 
Fig. 15 is very smooth and we think that the observed normal 
component of 12y is not produced by time-dependent fluctuations. 

It is concluded from the discussion in the Appendix that in 
crossing #4 there is a component of the magnetic field normal 
to the boundary; yet this crossing does not correspond to a 
rotational discontinuity. 

Another aim of the study of crossings #3 and #4 was to 
check the tailward drift hypothesis. For this purpose we 
consider these two crossings together. It is obvious from 
Fig. A^ and A^ that a change in the orientation of the boundary 
occurs between the crossings #3 and #4. For crossing #3 the 
average normal can be assumed to be roughly parallel to the 
X GSM axis (see the discussion in the Appendix) ; for crossing 
#4, the direction of the arrow is chosen. Fig. 16 summarizes 

our feeling about the shape of the boundary at the time of 
crossings #3 and #4, when we neglect the short period oscillations. 
The angles of 30° indicated in Fig. 16 are not claimed to be 
meaningful to better than +10°. The important point we want to 
emphasize is that the relative orientation of the two normals 
associated with the fact that #3 is observed before #4, implies 
that the structure AMB moves toward the tail. 
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What are the dimensions involved in this structure? First, 

we make the assumption that its tailward velocity is the 

velocity of the plasma flow (200 km/s) parallel to the upper- 

turbed boundary. The distance AB represents 2 minutes, or 

3 

120 x 200/6.4 10 'v 3.5 Re which can be compared with the 0.3 
Re obtained in the case of spike #2 (distance be in Fig. 10). 

To evaluate the thickness of the boundary we have to correct 
for the oscillations occurring during the crossing. The 
technique used is simple: from the slope of the unperturbed 
part of the crossing we evaluate graphically the duration of 
this crossing in the absence of oscillation. From Figs. 

A^ and A^ this gives about 2 sec for crossings #3 and #4. The 
angle between the velocity (200 km/s) and the boundary is +30° 
for the two crossings; we obtain a thickness of 200 km equivalent 
to 3 proton Larmor fadii. Alternatively, if this boundary 
thickness is considered to be a reasonable one, our result checks 
that the order of magnitude assumed for the velocity of the 
boundary was correct. So for the short or long period 
oscillations of the boundary (spike #2, or the pattern defined 
by crossings #3 and # 4), we reach identical conclusions: these 
oscillations are consistent with ripples on the boundary, drifting 
tailward with a velocity of the same order of magnitude as the 
plasma flow velocity. 
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5.3 Third Zone 

For the third zone of crossings after 1840, we tried to 
analyse the spikes labelled 5 and 6 in Fig. 4 by the technique 
already used for the spike #2 but this attempt was unsuccessful. 
Although no conclusion can be shown for this third zone we think 
that it can be useful to present a short account of the 
difficulties involved. Spike #5 appears very similar to spike #2 
even at a resolution of 7 sps. Its duration is about the same, 
namely 10 to 15 sec, but the eigenvectors computed from the 56 sps 
data for five different intervals were so scattered, with com- 
ponents normal to the boundary ranging from -15y to +20y, that 
we were unable to devise a scheme equivalent to that of Fig . ' 9 
for spike #5. For spike #6, the high resolution data (56 sps) 
presented in Fig. 17 shows that a study using the eigenvalue 
technique is meaningless: the fluctuations are too much rapid. 
Figure 17 represents the total field and the three components 
of the field in the satellite reference system. 

Two remarks about this spike can be useful. Till 1846:30, 

0G0-5 was in the magnetosheath and some seconds later the boundary 
swept partially over it. The first thing recorded by the satellite 
just outside the boundary was a very high magnetic field peaking 
at 70y during less than 2 sec at about 1846:33; such large 
amplitude fields were never recorded before 1800. After this 
time we have two more examples (see 1859 and 1913 in Fig. 4). 
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Another characteristic of the boundary at this time was 
the very high level of oscillations with periods less than 
1 sec. In this regard, it is interesting to compare Fig. 17 
with Figs. 7, 9, and A 2 which show crossings #1, #2, #3 and 
#4 with the same sampling rate. The main conclusion to be 
drawn from these very brief comments about the third zone is 
that owing to the great variability of the boundary at this time, 
one cannot assume this boundary to keep even a roughly constant 
orientation during the time required for one crossing, and so 
the usual analysis technique fails. 

6. DISCUSSION AND CONCLUSION 

We now summarize and discuss the main new findings of 
this study. ' 

Motion of the Boundary . From the observation of crossings 
we measured periods of several minutes and periods of less than 
one minute, even as low as 10 seconds. The two ranges do not 
seem to overlap. This is confirmed by the particle data. The 
long period oscillations are well within the period range 
reported by Hyde (1967), Anderson, et. al., (1968), Kaufmann 
and Konradi (1968), Smith and Davis (1970). They are shorter 
than those reported by Holzer, et. al., (1966). The Kelvin- 
Helmholtz instability is generally proposed as an explanation' 
of these oscillations (Southwood, 1968). We have shown on two 
examples (crossings #3 and #4) that the waveform associated with 
these oscillations drifted tailward with a velocity of the same 
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order of magnitude as the plasma flow velocity as would be 
required for a wave generated by the Kelvin-Helmholtz instability. 

The short period oscillations have not yet been reported in 
the literature, the main reason being that very often the 
sampling rate of the data did not permit the study of short 
periods. This was the case for Explorer 12 for instance (Hyde, 

1967, Kaufman and Konradi, 1968). In the case of IMP 2 the 

particle flux measurements were made every 40 sec; the fact 

that the spectrum of oscillation period obtained from IMP 2 

dropped sharply below 3 minutes (Anderson, et. al., 1968) 

confirm that the two ranges of oscillation namely several 

minutes and less than one minute are separated. We do not know 

if this implies two different generation mechanism, but we have 

shown that the 10 sec period ripple on the boundary responsible for 

the spike § 2 was also drifting tailward with a velocity similar to the 

plasma flow velocity. It would be interesting to know if the 

fluctuation of the. order of 10 seconds as represented by spikes 

#2, #5, and #6 have a different origin from the ones with less 

than 5 sec periods which appear in crossings #1 (Fig. 8), #2 

(Fig. 11), #3 (Fig. A l ) #4 (Fig. A 2 ) , and #6 (Fig. 17)., Let 

us recall that the proton cyclotron period in a 30y magnetic 

field is about 2 sec. The amplitude of the long period 

oscillations is generally much larger than the amplitude of the 

short period ones, except during the last part of the crossing 
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(between 1840 and 1906) . It should be emphasized that -this 
last sequence of crossings is characterized by: 

- the most southward orientation of the magnetosheath field 

- the most inward position of the boundary 

- the most fluctuating magnetic fields and particle .fluxes 
This is in contrast to the nearly steady boundary recorded at 
1700 where the magnetosheath field was northward. The oscil- 
lations as recorded by the particle experiment confirmed the 
qualitative description obtained from the analysis of the 
magnetic data. 

We assumed that the multiple crossings observed at the 
satellite are produced by the drifting toward the tail, of 
ripples on the boundary, with a velocity of the same order of 
magnitude as the average flow velocity of the plasma, (about 
200 km/sec) . The boundary thickness computed from this model 
was of 2 to 3 protons Larmor radii; if this thickness is 
considered to be a reasonable one, this is a rough check of the 
validity of our assumption on the velocity. The relative 
orientation of the boundary during succeeding crossings confirms 
that the motion was indeed tailward (a similar conclusion was 
reached by Kaufmann and Konradi, 1968). 

Microstructure of the Boundary 

We analyzed more than twenty crossings, four representative 
ones have been presented. Various signatures have been found 
as can be checked from Figs. 7 to 15. Moreover the evolution, 
as seen by the satellite, from one shape to the other is very 
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rapid: for example, there is a delay of only two minutes between 
the structures observed in Figs. 13 and 15. However, a fact 
remains constant through the whole interval: the difference 
vector of the magnetic fields on the two sides of the boundary 
is mainly vertical. The measurements being made at the equator, 
the difference vector for plasma flow velocity should be mainly 
horizontal. Therefore, we claim that the results cannot fit 
with a rotational discontinuity model of the boundary. This 
is confirmed by the fact that the amplitudes of the field on 
the two sides of the boundary always differ by at least 20%. 

On the other hand it has been shown also in section 5 that at 
some of the observed crossings^ normal components of the magnetic 
field are present and so these crossings cannot be interpreted 
as tangential discontinuities either. 

Normal Component 

The existence of a normal component of the magnetic field 
is probably the most important part of any study of the magneto- 
pause because it is customarily related to the open and closed 
models of the magnetosphere (Dungey 1961, Beard 1964). From 
a theoretical point of view, it has been shown by Siscoe (1966) 

that a uniform normal component over the whole magnetopause 

- 2 

should be roughly 10 y . Many previous attempts to determine the 
normal component of the magnetic field experimentally failed 
to obtain definitive results. One reason, as discussed in the 
Appendix, is that only high sampling rates permit the study 
of a highly time dependent magnetopause. However, in our opinion. 
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the main reason is that such a "definitive result" does not 
exist. First, our observations show that from one crossing to 
another very varying average normal components appear. In the 
four examples presented, these average normal components were 
mainly positive, hut we observed also negative or inward normal 
components. Secondly, inside each crossing the normal component 
could fluctuate by several gammas which is far from negligible 
when compared with a total field occasionally as low as lOy in 
the middle of the current sheet. Owing to -this variability, 
we cannot compare our results with the uniform normal component 
assumed by Siscoe (1966). These various normal components, 
which were never associated with a structure of rotational 
discontinuity, are undoubtedly related to the very high drag 
observed at this time and already reported in ARK. 

This study concerns a short interval of time and consequently 
does not pretend to give any statistical information. However, 
we had the opportunity to study a very perturbed state of the 
magnetopause and it is worthwhile to summarize the observations 
presented in ARK and in this paper. On March 27, 1968, following 
a reversal of the interplanetary magnetic field from northward 
to southward, some large amplitude oscillatory motions developed 
on the magnetosphere boundary; at the same time the very structure 
of this boundary was very variable and this variability 
increased with the southward component of the interplanetary 
field. Associated with these phenomena, variable components 
of the magnetic field normal to the boundary were recorded and 
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produced a large transfer of momentum across the magnetopause 
leading to the erosion of the dayside magnetosphere. It must 
be emphasized that during all this period of time, the observed 
local characteristics of the boundary, including the amplitude 
and orientation of the normal component as well as the counting 
rate of energetic electrons on either side of the boundary, were 
extremely variable and thus did not fit with any steady open or 
closed model of the magnetosphere. 
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APPENDIX 


For the detailed study of individual crossings, the first 
problem is to display the variations of the magnetic field in a 
convenient reference system. In particular we wish to choose 
the x and y directions to be in the plane of the boundary and 
z normal to the boundary. The technique for determining this 
reference system is now very well known (Sonnerup and Cahill 1967) 
Since B must satisfy 

V ' B = 0, ( A . 1) 

if we assume a plane homogeneous boundary such that 

9B 9B 

3 = 9 = 0 (A. 2) 

x y 

then 


9 B 

- 0 (A. 3) 

9 

z 

2 

If Cf is the variance of the projection, on a direction _n, of 

l 

a sample of N measurements of the magnetic field, 
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£ 

i=l 




(A. 4) 
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2 - 

then Eq. (A. 3) implies that O should be zero when is parallel 

to the normal to the boundary, but in practice this will not 

be entirely true due to the presence of waves and because the 

boundary is not planar nor homogeneous. ■ The maximum and minimum 
2 

values of O are obtained by computing the eigenvalues of the 
matrix 


M 


a3 


= B B. 
a 


- B B. 

a ( 


where B and B c are components of an individual measured field 

a 3 

vector in the reference system of the satellite. 

The vector associated with the smallest eigenvalue is 
normal to the boundary. The vectors jc. ahd y associated respectively 
with the largest and intermediate eigenvalues, are in- the 
plane of the boundary. If the current inside the boundary is 
unidirectional, 2 : is perpendicular to the current vector and 
then x anc * .£ are undetermined. 

Some remarks should be made: the method is based on the 
validity of Eqs. (A. 2) and (A. 3) which requires a plane and steady 
state boundary during the time of the analysis. A crossing as 
seen by 0G0 5 lasts about 10 seconds and we have shown evidence 
for important variations of the orientation of the boundary 
within such a time interval. Thus, at least in our application, 

Eqs. (A. 2) and (A. 3) are not generally satisfied over an entire 
crossing. It is, therefore, important in applying the eigenvalue 
technique to boundary crossings to select for analysis, intervals 
short enough that the boundary remains relatively steady. In 
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this appendix we illustrate for two crossings the effect of 
changing the interval of analysis across a boundary traversal 
and of changing the sampling rate of the data used in the analysis. 
For the UCLA QGO-5 fluxgate, changes in sampling rate are 
equivalent to averaging the data points over different sample 
periods because the data is filtered to obey the Nyquist 
criterion at each sampling rate. 

First, let us make an eigenvalue analysis for crossing #3. 

An enlarged plot of the B^ component for this crossing measured 
at a rate of 56 sps is shown in Fig. A1 . Eigenvalue computations 
have been processed for eleven time intervals shown on this 
figure. For the intervals labelled to a rate of 56 sps 

was used; for the intervals labelled n^ and n.^, 1 sps was used. 

Table 1 lists the following information for these eigenvalue 
computation's: the number of samples used to perform the eigen- 
value computation, the average component B z of the magnetic field 

2 

normal to the boundary as .well as the variance a (Eq. (A4))of 
the normal component of the field around this average value; 

(this of course equals the smallest eigenvalue) . The ratio of 
the largest eigenvalue to the intermediate one as well as the 
ratio of the intermediate to the smallest one are also shown 
for each interval. 

Table 1 shows several facts about the eigenvalue analysis. 
First, the field along the computed normal is quite variable 
ranging from 10 to -20y. Secondly, the variance of the field 
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2 

along the normal can be quite small, as low as . 03y and is 
largest for the intervals covering the longest period of time. 
Thirdly, the three eigenvalues generally have quite distinct 
magnitudes. We shall return to the first point later. The 
second point implies that the boundary orientation does not 
change much over short periods but can move significantly over 
long periods. The third point means that the direction of' 
maximum field variation (the x axis) and of minimum field 
variation (the boundary normal) are both well defined. However, 
the 56 sps data generally give a better separation of the 
intermediate and minimum eigenvalues than does 1 sps data, and 
thus we expect the 56 sps data defines the boundary normal better. 

The directions of the normal to the boundary (arrows), and 
of the x axis (circle) obtained in these eleven computations are 
shown at the upper left of Fig. A1 . The representation is the 
same as in Fig. 6, namely each point (arrow or circle) represents 
the projection in the GSM equatorial plane of the tip of a 
unitary eigenvector. Circles corresponding to angles of 10°, 

30° and 90° with the GSM Z axis are shown. The dashed arrows 
(N^, Ng , and Nj.) represent vectors below the XY GSM plane. 
Except for these four vectors, all the eigenvectors in Fig. A1 
have a positive component along the GSM Z axis. 

The variation of the boundary orientation during this 
crossing can be studied from the evolution of the z_ eigenvectors 
(arrows). These vectors have been grouped in three sets (shaded 
areas) : and correspond to the first 1 second part of the 
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crossing around 1817:08 UT ; N, , N Q and correspond to the last 
part, 2 seconds around 1817:12; the intermediate part lasting 
3 seconds is represented by N^, N,. and Ng. The £ vector is 
computed from the whole interval and so represents some time 
average of the boundary normal direction. However, one must not 
forget that the time spent by the satellite inside a region 
of the boundary at a particular time depends on the angle between 
the instantaneous normal and the velocity vector of this part 
of the boundary (a velocity vector parallel to the normal gives 
the shortest crossing, whereas a velocity vector parallel to 
the boundary gives the largest one) . Consequently the time- 
averaged normal (as the arrow) can be different from the 
space-averaged normal we are interested in. For this particular 
crossing (#3) the latter normal is probably not far from the 
X axis. The results obtained from the 1 sps data are rather 
different from the results obtained from the 56 sps data even 
when the interval is exactly the same (n.^ and for instance) 
and we show them only to illustrate the advantage of the high 
time resolution data. 

The largest eigenvectors (circles) are on the average not 
very far from the Z GSM axis (see and n.^) a result 

which would be obtained from a current sheet flowing on the 
average very roughly parallel to the GSM equatorial plane. 

However, the structure of this current sheet is rather complicated 
as can be checked from the scattering of the largest eigenvectors 
obtained from short time samples (^ to N^). 
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Let us now examine the eigenvalue analysis for the crossing 
#4. An enlarged plot of the variation of the component 
measured at a rate of 36 sps is shown in Fig. A2 . Ten intervals 
chosen for the eigenvalue computation are delineated and the 
directions of the corresponding axes jc and _z are shown in the 
upper right of the figure. 

Table 2 presents the information about this set of eigen- 
value computations in the same format as used in Table 1. As 
before the component of the field normal to the boundary varies 
over the different analysis intervals, ranging from 4.4 to 13. 6y. 
Again the ratio of the intermediate to the smallest eigenvalue 
illustrates that the 1 sps data can only poorly define the 
direction of the boundary normal. Also, as before the longer 
intervals are associated with larger variances. 

Examining Fig. A2 , we see that the eigenvectors are less 
scattered than those in crossing #3. In fact the main scatter 
is due to ng , n^ and Nj • The calculation of ng and n^ used 
1 sps data which ^as discussed above^poorly define the normal, 

-and the interval over which was calculated includes a bump 
which appears to be due to a short period fluctuation with a 
period of the order of seconds of the same type which produced 
an intermediate zone in crossing #3. If these three normals^are 
ignored, the grouping of the calculated normals around some 
average normal direction is apparent. 

The computed x-axes (circles) lie as expected, roughly in 
a plane perpendicular to this average normal. (The intersection 



of this plane with the unit sphere is indicated by the dashed 
line.) Here, again, the x-axes corresponding to the long 
duration samples , n g , n^ and n 1Q are closest to the Z GSM 
axis which is consistent with a current sheet flowing approxi- 
mately in the GSM equatorial plane. The scatter of the x-vectors 
computed over the shorter Intervals reflects the complexity 
of the detailed structure of the current sheet. 

From an examination of these two crossings we see first 
that we generally cannot determine consistent boundary normals 
with the 1 sps data. Furthermore, even the 56 sps data may 
give misleading normals at times when the boundary oscillates 
during an analysis interval. Thus we must use the results of 
the eigenvalue technique with caution^ requiring a consistent 
pattern of normals over short periods within the boundary if 
we are to obtain meaningful results. High data rates such as 
available on 0G0-5 are essential for this purpose, although 
low sampling rates certainly are adequate to give a general 
picture of the field variation across the boundary (as in Figs. 

8, 10, 12 and 14) . 

With this in mind we can interpret the extent to which the 

field component normal to the boundary as derived in Tables 1 

and 2 is meaningful. First, we should ignore all the components 

derived from the 1 sps data. In the crossing #3 the various 

B components (Table 1) are very variable and both positive and 
2 

negative so we cannot claim a steady normal component exists 
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during this crossing. On the other haiid , if we discard the 
1 sps data, the normal components obtained for the various time 
intervals in the crossing #4 (Table 2) are rather constant of 
the order of 10y* The only exception is the interval N 2 which 
is a perturbed part of the crossing. The rather steady normal 
component during the remainder of the crossing shows that this 
measured normal component is not a consequence of fluctuations 
in the orientation of the boundary. As the error in zero-level 
discussed in section 2 is less than 2y for each axis, a zero- 
level error cannot explain such a lOy normal component either. 
Thus, we see that at times there can exist large field- components 
normal to the magnetopause. 
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Captions 


Fig . X . 


Fig. 2. 


Fig . 3 . 


Fig. 4. 


Fig. 5. 


Fig. 6. 


Projection of the orbit of 0G0 5 in the GSM 
equatorial plane at the time of the observations. 

The extreme extrapolated • positions of -the. .boundary 
are shown. 

Basic relations for tangential and rotational 
discontinuities. The projection of the magnetic 
field vectors on the two sides of the discontinuity, 
in the plane of this discontinuity are shown on 
the right. 

Variation of the magnetic field versus universal 
time an March 27, 1968. The GSM reference system 
is used; B^, refers to the total field. 

Variation of the magnetic field and particle 
flux versus time from 1719 to 1919 UT on March 27, 
1968. The numbers from 1 to 6 refer to crossings 
which will be studied in Section 5. The particle 

A 2 

flux unit is 10 H particles/cm /sec/ster. See 
text for details. 

Schematic models of an equatorial section of the 
boundary . 

Projection on the X Y GSM plane of the tips of 
unitary vectors parallel to the largest eigen- 
vectors (x in the text, open circles in the drawing) 
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and the smallest eigenvectors (_z in the text* 
arrows in the drawing) for several crossings. 

2 ^ and x^ refer to the first crossing studied 
(see text) . z and x refer to the unperturbed 
boundary." The circles indicate the locus of 
constant angle 30°, 45° and 90° with the Z, GSM 
axis . 

Fig. 7. Variation of the magnetic field through crossing 

as seen in the reference system xyz of the eigen- 
vectors. The sampling rate is 56 sps. 

Fig. 8. Evolution of the tip of the projection of the 

magnetic field in the plane xy of the boundary 
for crossing #1. Points are 1 sec apart. 

Fig. 9 Variation of the magnetic field during spike #2 

as shown in the reference system of the satellite. 

is the total field. A aero base line is drawn 
for each component. A sampling rate of 56 sps 
is used. 

Fig. 10. Shape of the boundary at the time of #2. The 

circled numbers are the values of B^, computed 
from 1 sps data; the associated short arrows give 
the direction of the field projection in the zy 
plane. On the left, are shown various relevant 
directions as seen in the GSM equatorial plane 
(see text) . 
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Fig. 11. 


Fig. 12. 


Fig. 13. 


Fig. 14. 


Fig. 15. 


Variation of the magnetic field through crossing 
b (Fig. 10) as seen in the reference system xyz 
of the eigenvectors. 

Variation of the magnetic field component in the 
eigenvector reference system xyz for interval n^ 
(crossing #3, Fig. Al) . Points are 1 sec apart. 
Evolution of the tip of the projection of the 
magnetic field in the plane xy of the boundary 

for interval (crossing #3, Fig. Al) . The 

_ 2 

points are 1.7 x 10 sec apart (56 sps) . The 
dashed lines are drawn to emphasize the regions 
where VB is parallel to IS. The average value of 

(3.5 10 y) for the interval-AB is less than the 
jiigitization error of the instrument and so is 
equivalent to zero. 

Variation of the magnetic field components in 

the eigenvector reference system for interval 

n Q (crossing // 4, Fig. A2) . Points are 1 sec 
o 

ap ar t . 

Evolution of the tip of the projection of the 

magnetic field in the plane xy of the boundary 

for interval (crossing #4, Fig. A2) . The 

-2 

points are 1.7 x 10 sec apart (56 sps). 
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Fig. 16. Sketch of the boundary at the time of the 

crossings #3 and #4. Not to scale. 

Fig. 17. Variation of the magnetic field through the 

spike #6 (Fig. 4). The satellite reference 
system 1, 2, 3 is used, as in Fig. 9 . A zero 
base line is drawn for each component. A 
sampling rate of 56 sps is used. 

Fig. Al. Variation of the component (reference system 

of the satellite) versus time, through crossing 
#3. The sampling rate is 56 sps. The intervals 
labelled to n^ were used for eigenvalue 
computations. The resulting directions of the 
largest (open circles) and smallest (arrows) 
eigenvectors are represented at the top left of 
the figure (same system of representation as in 
Fig. 6 ). x and z refer to the unperturbed 
boundary . 

Fig. A2. Variations of the B^ component (reference system 

of the satellite) versus time through crossing #4. 
The sampling rate is 56 sps. The intervals 
labelled to n^Q were used for eigenvalue 
computations. The resulting directions of the 
largest (open circles) and smallest (arrows) 
eigenvectors are represented at the top right 
of the figure (same system of representation as 
in'Fig. 6 ). x and z refer to the unperturbed 
boundary. 
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Table 1,. 


Table 2. 


Results of the eigenvalue computation for several 
time intervals during the crossing #3. 

Results of the eigenvalue computation for several 
time intervals during the crossing #4. 
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Table 1 


Eigenvalue Ratios 



Number 


o 

Maximum 

Intermediate 

Interval 

of' Samples 

b z (y) 

Z 

a 

Int ermediate 

Minimum 

N 1 

625 

+ 2.6 

7 . 2 

18.8 

5.9 

N 2 

73 

-10.2 

0.1 

9.1 

24.7 

N 3 

36 

-6,1 

0.03 

56.9 

17.1 

N 4 

58 

-22.4 

0.1 

6.8 

2.9 

N S 

51 

-10.0 

0.04 

12.9 

5.7 

N 6 

51 

+7.3 

0.2 

6.5 

4.0 

N 7 

100 

-0.5 

1.5 

6.3 

33.6 

N 8 

50 

-0.7 

2.0 

7.3 

16.7 

N 9 

50 

+0 . 4 

0.3 

5.4 

8.8 

n io 

8 

+10 . 3 

7.7 

31.6 

2.8 

n ll 

11 

+4.7 

5.9 

29.7 

4 . 4 
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Table 2 






Eig envalue 

Ratios 


Number 



Maximum 

Intermedi 

Interval 

of Samples 

b z (y) 

2 

a 

Intermediate 

Minimum 

N 1 

109 

13.6 

0.1 

11.2 

198.4 

*2 

45 

4.4 

0.1 

5.1 

19.1 

N 3 

67 

9 . 2 

0.8 

12.8 

11.4 

N 4 

221 

11.1 

1.9 

7.3 

43.4 

N 5 

36 

10.8 

0.03 

135.2 

13.5 

N 6 

59 

12.4 

0.1 

38.8 

21.8 

N 7 

669 

8.8 

3.7 

23.4 

15.3 

n 8 

12 

13.2 

6.0 

90.8 

2.3 


14 

12.2 

9.7 

104.7 

1.3 


n 10 


7 


14.5 


4.3 


56.7 


3.9 
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